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Heat Transfer: Exercises Preface 


Preface 


Worked examples are a necessary element to any textbook in the sciences, because they reinforce the 
theory (i.e. the principles, concepts and methods). Once the theory has been understood, well chosen 
examples can be used, with modification, as a template to solve more complex, or similar problems. 


This work book contains examples and full solutions to go with the text of our e-book (Heat Transfer, 
by Long and Sayma). The subject matter corresponds to the five chapters of our book: Introduction to 
Heat Transfer, Conduction, Convection, Heat Exchangers and Radiation. They have been carefully 
chosen with the above statement in mind. Whilst compiling these examples we were very much aware 
of the need to make them relevant to mechanical engineering students. Consequently many of the 
problems are taken from questions that have or may arise in a typical design process. The level of 
difficulty ranges from the very simple to challenging. Where appropriate, comments have been added 
which will hopefully allow the reader to occasionally learn something extra. We hope you benefit 
from following the solutions and would welcome your comments. 


Christopher Long 
Naser Sayma 
Brighton, UK, February 2010 
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1. Introduction 


Example 1.1 


The wall of a house, 7 m wide and 6 m high is made from 0.3 m thick brick with k =0.6W/mkK . 


The surface temperature on the inside of the wall is 16°C and that on the outside is 6°C. Find the heat 
flux through the wall and the total heat loss through it. 


inside 


outside 


T,= 6°C T, = 16°C 


Solution: 


For one-dimensional steady state conduction: 


dT ok 
af 2-8 =f 
q 7, ) 
0.6 
ee gb) D0 ti 
q met ) m 


O = qA=-20x(6x7)=-840W 


The minus sign indicates heat flux from inside to outside. 
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Example 1.2 


A 20 mm diameter copper pipe is used to carry heated water, the external surface of the pipe is 
subjected to a convective heat transfer coefficient of h = 6 W/m°K , find the heat loss by convection 


per metre length of the pipe when the external surface temperature is 80°C and the surroundings are at 
20°C. Assuming black body radiation what is the heat loss by radiation? 


T,= 20°C 


T, = 80°C 


20 mm diameter 


Solution 

om = A(T, ~T, )= 680-20) = 360 W/m? 

For | metre length of the pipe: 

Ono = Veonr = Vecomy X 244 = 360x 2x 2x 0.01 = 22.6W/m 
For radiation, assuming black body behaviour: 

Gua = (0, -T?) 

Qa = 5-67 10° (3534 — 293°) 

Drag = 462W | m? 

For | metre length of the pipe 

Ont = Ung A = 462 x 2x 1x 0.01 = 29.1 W/m? 


A value of h = 6 W/m’ K is representative of free convection from a tube of this diameter. The heat 
loss by (black-body) radiation is seen to be comparable to that by convection. 
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Example 1.3 


A plate 0.3 m long and 0.1 m wide, with a thickness of 12 mm is made from stainless steel ( 

k =16W/mkK ), the top surface is exposed to an airstream of temperature 20°C. In an experiment, 
the plate is heated by an electrical heater (also 0.3 m by 0.1 m) positioned on the underside of the plate and 
the temperature of the plate adjacent to the heater is maintained at 100°C. A voltmeter and ammeter are 
connected to the heater and these read 200 V and 0.25 A, respectively. Assuming that the plate is 
perfectly insulated on all sides except the top surface, what is the convective heat transfer coefficient? 


W=0.1m 


k=16W/mK 


t=0.012m 


insulated -» “\ NN NN 


Solution 
Heat flux equals power supplied to electric heater divided by the exposed surface area: 


_VxI_ Vx _ 200x0.25 


= = = 1666.7 W/m? 
A WxL 0.1x 0.3 


This will equal the conducted heat through the plate: 


k 
q=—(T -T,) 


7 =7, 4 = 199 1966-7 0-012) _ 99 a50¢ (371.75 K) 
tk 16 


The conducted heat will be transferred by convection and radiation at the surface: 
q= HT, -1,)+ olf? -7/) 


q-o(T;'-T;) 1666.7 -5.67x10°*(371.75* - 293°) 
Ca 371.75 — 293 


=12.7W/m°K 
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Example 1.4 


An electronic component dissipates 0.38 Watts through a heat sink by convection and radiation (black 
body) into surrounds at 20°C. What is the surface temperature of the heat sink if the convective heat 
transfer coefficient is 6 W/m’ K, and the heat sink has an effective area of 0.001 m’? 


Solution 

g=£=i(r, -7,)+o(t' 72) 

038 _ 6(r, —293)+5.67x107(T* — 293°) 
0.001 


5.67x10°T? +67, — 2555.9 =0 
This equation needs to be solved numerically. Newton-Raphson’s method will be used here: 


f =5.67x10°T* +67, — 2555.9 


W opebs. 10*T? +6 
aT. 


Ss 


5.67x10°T? +67, — 2555.9 


° df ° 22.687’ +6 
dT, 


ha = T" f = T"” 


Start iterations with T° =300 K 


-8 4 
T! =300 5.67 x10 x300° + 6x 300 2555.9 _ 35446 K 


22.68 x 300° +6 


-8 4 _ 
T? = 324.46 5.67x10~° x 324.46" +6x 324.46 2555.9 _ 353K 


22.68 x 324.46° +6 
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Introduction 


The difference between the last two iterations is small, so: 


T° =323 K =50°C 


The value of 300 K as a temperature to begin the iteration has no particular significance other than 
being above the ambient temperature. 
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2. Conduction 


Example 2.1 


Using an appropriate control volume show that the time dependent conduction equation in cylindrical 
coordinates for a material with constant thermal conductivity, density and specific heat is given by: 


OT ier 7. er 
+ + = 
or? ror Gz” act 


Were @ = s is the thermal diffusivity. 
pc 


Solution 


0.25 z 


2,5 


Consider a heat balance on an annular control volume as shown the figure above. The heat balance in 
the control volume is given by: 


Heat in + Heat out = rate of change of internal energy 


Ou 
QO, + QO, ae er =U a oe (2.1) 
t 
Og, = Q, +25, 
or 
0.45: = QO. + 2 5, 
Oz 
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u=mcl 
Substituting in equation 2.1: 


eo ee ee (ee 
or Oz Ot 


Fourier’s law in the normal direction of the outward normal n: 


CG 

A On 

eT Ee Er ae (A=27réz) 
Or Or 

6-387 24 (A=2z7roér) 
Oz Oz 


Equation 2.1 becomes 


loz = MC a (2.3) 


a ee ae eee rr a 
0 Or 0 0 Ot 


ia Z Z 
Noting that the mass of the control volume is given by: 


m=p2rrdérdz Equation 2.3 becomes 


Shir Zhore Zhe a = pepe 
or or Oz Oz Ot 


Dividing by r, noting that r can be taken outside the brackets in the second term because it is not a 
function of z. Also dividing by k since the thermal conductivity is constant: 


10 a OT pcoT 
r +f = 
ror| or Oz’ k Ot 


k : 
Using the definition of the thermal diffusivity: @ =— and expanding the first term using the product rule: 
pc 


which gives the required outcome: 


2 oT To Or _ 158 


r—+ + = 
ror| Or? Or or z> a Ot 
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oT er .@T  1oF 
+ + = 
or? ror 62” act 


Example 2.2 


An industrial freezer is designed to operate with an internal air temperature of -20°C when the external 
air temperature is 25°C and the internal and external heat transfer coefficients are 12 W/m’ K and 8 
W/m’ K, respectively. The walls of the freezer are composite construction, comprising of an inner 
layer of plastic (k = 1 W/m K, and thickness of 3 mm), and an outer layer of stainless steel (k = 16 
W/m K, and thickness of 1 mm). Sandwiched between these two layers is a layer of insulation material 
with k = 0.07 W/m K. Find the width of the insulation that is required to reduce the convective heat 
loss to 15 W/m”. 
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insulation k; = 0.07 W/m K 


ease outside 
ee T= 25°C 
— 2 
hi= 12 Wim? K \ a 
N q 
3mm plastic 
k,= 1Wi/m K 1mm Stainless 
steel, k, = 16W/mK 
Solution 


q=UAT where U is the overall heat transfer coefficient given by: 


a 15 
AT 25-(-20) 


= 0.333W/m?K 


L 
p=<g1— L oF, 1 Wool Piece] 
0.333 |h, k, k, Ah, 0.333 [12 1 16 8 


L, =0.195m (195 mm) 
Example 2.3 


Water flows through a cast steel pipe (k = 50 W/m K) with an outer diameter of 104 mm and 2 mm 
wall thickness. 
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i. Calculate the heat loss by convection and conduction per metre length of uninsulated pipe 
when the water temperature is 15°C, the outside air temperature is -10°C, the water side heat 
transfer coefficient is 30 kW/m” K and the outside heat transfer coefficient is 20 W/m’ K. 

ii. Calculate the corresponding heat loss when the pipe is lagged with insulation having an 
outer diameter of 300 mm, and thermal conductivity of k = 0.05 W/m K. 


Solution 
T, 
T, = -10°C 
h, = 20 W /m°K 
T, =15°C 
h, = 30 kW /m°*K 
k=50W /mK 
Plain pipe 
Q 
=2arLh{T. -T, > T-T = 
QO 1 af i :) i 1 on r,Lh, 
QaLk\T, —T. 
O= ( 1 ») > i, 7 = QO 
In(r, /7,) 2a Lk /\n(r,/7,) 
Q 
=2azr,Lh |T, -T > T, -T.= 
QO 2 Aik 2 ’) 2 (a) In r,Lh, 


Adding the three equations on the right column which eliminates the wall temperatures gives: 


2aL(T, -T,) 
o= 1 In(r. /r,) 1 

+2 ty 
hr, k hr, 


oO 


Q 2n(15-(-10)) 
- =163.3W/ 
L I, In(0.052/0.05) : 
30000 x 0.05 50 20x 0.052 
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Insulated pipe 


O_ 2n(T, -T,) 
L 1 In(x,/r) rn In(@y,/r,) 1 
hr, k Bins hr, 
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2n(15—(-10)) 
1 r In(0.052/0.05) , In0.15/0.052) 1 
30000 x 0.05 50 0.05 20x 0.15 


Q =7.3W/m 
L 


For the plain pipe, the heat loss is governed by the convective heat transfer coefficient on the outside, 
which provides the highest thermal resistance. For the insulated pipe, the insulation provides the 
higher thermal resistance and this layer governs the overall heat loss. 


Example 2.4 


Water at 80°C is pumped through 100 m of stainless steel pipe, k = 16 W/m K of inner and outer radii 
47 mm and 50 mm respectively. The heat transfer coefficient due to water is 2000 W/m’ K. The outer 
surface of the pipe loses heat by convection to air at 20°C and the heat transfer coefficient is 200 W/m” 
K. Calculate the heat flow through the pipe. Also calculate the heat flow through the pipe when a layer 
of insulation, k = 0.1 W/m K and 50 mm radial thickness is wrapped around the pipe. 


Solution 
The equation for heat flow through a pipe per unit length was developed in Example 2.3: 


2aL(T, -T,) 
1 n In(r, in), 1 
hin, k hr, 


L oO 


O= 


Hence substituting into this equation: 


2m x 100(80 — 20) 
1 A In(50/47) st 
0.047 x 2000 16 0.05 x 200 


O= = 0.329 x 10°W 


For the case with insulation, we also use the equation from Example 2.3 


2aL(T, -T,) 
o- 7 tain) in@eiey. 1 
aw 2 1 4 3 2 4 
hin, k he h,r; 
O= 2 x 100(80 — 20) 2 eabeaa 


1 . In(50/47) | In(100/50) I 
0.047 x 2000 16 0.1 0.1x 200 


Notice that with insulation, the thermal resistance of the insulator dominates the heat flow, so in the equation 
above, if we retain the thermal resistance for the insulation and ignore all the other terms, we obtain: 
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O- 2aL(T,-T,) _2xx100(80-20) 
~ In(r,/r,) ~———«sIn(100/50) 
k 0.1 


ins 


=5.44x10°W 


This has less than 1% error compared with the full thermal resistance. 
Example 2.5 


A diagram of a heat sink to be used in an electronic application is shown below. There are a total of 9 
aluminium fins (k = 175 W/m K, C = 900 J/kg K, p = 2700kg/m’) of rectangular cross-section, 
each 60 mm long, 40 mm wide and | mm thick. The spacing between adjacent fins, s, is 3 mm. The 
temperature of the base of the heat sink has a maximum design value of T, = 60°C , when the external 
air temperature 7’, is 20°C. Under these conditions, the external heat transfer coefficient h is 12 W/m” 


K. The fin may be assumed to be sufficiently thin so that the heat transfer from the tip can be 
neglected. The surface temperature T, at a distance, x, from the base of the fin is given by: 


T, —T, )coshm(L — x 
T-T, = ( pad } ) where m* = coe and A_ is the cross sectional area. 
sinh mL A 


(og 


Determine the total convective heat transfer from the heat sink, the fin effectiveness and the 
fin efficiency. 


7,=2rc 
h=12W/m?K 
E 
= 
oO 
wo 
=, 
T, = 60°C 
N = 9 fins 
K=175WimK 
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Conduction 


Solution 


Total heat fluxed is that from the un-finned surface plus the heat flux from the fins. 
QO=9, +0, 

Q, = A,h(T, -T,) = wxs(N-D)a(T, -7,) 

Q,, = 0.04 x 0.003(9 —1))x 12(60 — 20) = 0.461 W 


For a single fin: 


aT 
i cies {¢ dx ia 


Where A. is the cross sectional area of each fin 


Since 
(7, -T, )cosh m(L — x) 


T-T, = - 
sinh mL 


iB \ 
3 It’s only an 


opportunity if 


you act on it 


IKEA.SE/STUDENT 


600Z 'A'g SweIsAs WIy] JEW] @ 
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Then 
dT = a ae T,) 
dx cosh mL 
Thus 
aT —sinhmL 
.=—kA | — =—-kA m\T, —T, 
2, {4 a if cosh mL ( _ ') 


OQ, =kA.m(T, —T, )tanh(mL) = (hpkA, (7, —T, )tanh(mL) 
Since 
hP \’ 
pee 
a) 
P =2(w+1) = 2(0.04 x 0.001) = 0.082 m 


A, =wxt=0.04x 0.0001 = 40x10 m? 


1 


ee * APBSE 7 
es) u 


-( 12x 0.082 
mL = 11.856 0.06 = 0.7113 


tanh(mL) = tanh(0.7113) = 0.6115 


O, =(120.082x175x 40x10)” x (60-20) 0.6115 = 2.03 W / fin 
So total heat flow: 

O=0, +0, =0.461+9x2.03 =18.7 W 

Finn effectiveness 


: Fin heat transfer rate _ QO, 
Heat transfer rate that would occur in the absence of the fin AA, (7 oe i) 


an 


— 2.03 
12 x40x10~6(60 — 20) 


= 106 
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Fin efficiency: 


Actual heat transfer through the fin 


= Heat that would be transferred if all the fin area were at the base temperature 
n = 
WAT, -T, 


A, =wLh+wlh+Lt+ Lt =2L(w+t) 
A, = 2x 0.06(0.04 + 0.001) = 4.92 x10° m’ 


2.03 


= = 0.86 
12x 4.92 x10°(60—20) 


7 jin 


Example 2.6 
For the fin of example 4.5, a fan was used to improve the thermal performance, and as a result, the 
heat transfer coefficient is increased to 40 W/m’ K. Justify the use of the lumped mass approximation 


to predict the rate of change of temperature with time. Using the lumped mass approximation given 
below, calculate the time taken, 7 , for the heat sink to cool from 60°C to 30°C. 


Solution 


Consider a single fin (the length scale L for the Biot number is half the thickness t/2) 


p _ hh _hxt/2 _ 40x0.0005 


: ~107* 
k k 175 


Since B, << 1, we can use he “lumped mass” model approximation. 


(r-T,) ea 
= exp— 
(7 -T,) mC 
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m= pAt/2 


T-T _ 
a On f _ _ 2700x900 x 0.001 30-20 AP acesuils 


aa 2x40 60 — 20 


Example 2.7 


The figure below shows part of a set of radial aluminium fins (k = 180 W/m K) that are to be fitted to 
a small air compressor. The device dissipates 1 kW by convecting to the surrounding air which is at 
20°C. Each fin is 100 mm long, 30 mm high and 5 mm thick. The tip of each fin may be assumed to be 
adiabatic and a heat transfer coefficient of h = 15 W/m’ K acts over the remaining surfaces. 


Estimate the number of fins required to ensure the base temperature does not exceed 120°C 
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Conduction 


Solution 
Consider a single fin: 


P =2(w+t) = 2(0..005 x 0.03) = 0.07 m 


A, = wxt =0.005 x 0.03 =150x10° m? 


(ae e 15x 0.07 -| = 6.2361 m* 
kA, 180x150x 107 


mL = 6.2361x 0.1 = 0.62361 


tanh(mL) = 0.5536 
OQ = (APkA, ae (7 io a f )tanh(mL) (From example 2.5) 


O, =(15x0.07x180x150x 10°)” x (120-20)x 0.5536 = 9.32 W 


So for 1 kW, the total number of fins required: 


N= 4000 108 
9,32 
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Example 2.8 


An air temperature probe may be analysed as a fin. Calculate the temperature recorded by a probe of length 
L=20 mm, k= 19 W/m K, D=3 mm, when there is an external heat transfer coefficient of h = 50 W/m’K, 
an actual air temperature of 50°C and the surface temperature at the base of the probe is 60°C. 


Solution 


The error should be zero whenT,,, = T,, . The temperature distribution along the length of the probe 


(from the full fin equation) is given by: 


hi. 
7p. coshm(L—x)+ a sinh m(L — x) 


h, 
is coshmL +—” sinh mL 
mk 
1/2 
m-( A=2zD’/4, P=xzD 


At the tip, x = L, the temperature is given by (cosh(0) =1, sinh(0) = 0): 


Be =. 1 
“4 a " =¢ 
6? coshmL + —? sinh mL 
mk 


Where ¢ is the dimensionless error:, 


= 0, t= 1, (no error) 
tp 
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¢=1, T, =T, (large error) 


For L=20mm, k=19W/mkK, D=3mm, h=h,, =50W/m’ K 


tip 


T,, =50°C, T, = 60°C 


A=2zD* /4, P=xzD 


1/2 1/2 1/2 1/2 
mee) eee] ie (=) z es) = 59.235 m” 
kA kz D kD 19 x 0.003 


mL = 59.235 x 0.02 = 1.185 


t = 2 = 0.0444 
mk 359.235%19 
5 ee a 1 


= 0.539 


T,-T, cosh(1.185)+0.0444 x sinh(1.185) 
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T, 


tip 


205307 =7 er. 


T,,, = 0.539(60—50)+ 50 = 55.39°C 


tip 
Hence error = 5.39°C 
Example 2.9 


A design of an apartment block at a ski resort requires a balcony projecting from each of the 350 
separate apartments. The walls of the building are 0.3 m wide and made from a material with k = 1 
W/m K. Use the fin approximation to examine the implications on the heat transfer for two separate 
suggestions for this design. In each case, the balcony projects 2 m from the building and has a length 
(parallel to the wall) of 4 m. Assume an inside temperature of 20°C and an outside temperature of - 
5°C; the overall (convective + radiative) heat transfer coefficient on the inside of the building is 8 
W/m’ K and on that on the outside of the building is 20 W/m? K 


a) A balcony constructed from solid concrete and built into the wall, 0.2 m thick, k = 2 W/m K. 
b) A balcony suspended from 3 steel beams, k = 40 W/m K, built into the wall and projecting out by 2 


m each of effective cross sectional area A, = 0.01 m, perimetre P = 0.6m (The actual floor of 


the balcony in this case may be considered to be insulated from the wall 
c) No balcony. 


Solution 


W=2m 


Wall 
k,=1W/mK 


balcony 
k,=2W/mK 


T, =-5°C 


h, =20W/m?K 


L=0.3m 


a) For the concrete balcony 


h,t/2 


oO 


k, 


Treat the solid balcony as a fin B, = 
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Not that Bi is not << 1, thus 2D analysis would be more accurate. However, treating it as a fin will 
give an acceptable result for the purpose of a quick calculation. 


P=2(H +t) =2(44+0.2)=8.4m 
A, =Hxt=4x0.2=0.8m’ 


To decide if the fin is infinite, we need to evaluate mL (which is in fact in the notation used here is mW) 


1/2 1/2 
mw (He) w (xs) x2=20.5 


KA 2x0.8 


This is large enough to justify the use of the fin infinite equation. 


O, = (h, Pk, A, eee -T,) 


1/2 
1 h, Pk 
t= Lh Phd), -17,)-[ MEE 1) 0) 


Cc Cc 


Also assuming 1-D conduction through the wall: 


q, =h, 7, -T) (2) 


k 
q) => -h) (3) 


Adding equations 1, 2 and 3 and rearranging: 


Cas) 
1/2 
1 L(A, 
+—+ 
h, k, [4 | 


This assumes 1D heat flow through the wall, the concrete balcony having a larger k than the wall may 


(4) 


qy = 


introduce some 2-D effects. 
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From (4) 


20 -(-5) 
qs = > =77.2 Wim 


1 03 0.8 
+——+ 
Ss. 2 20x8.4x2 


Compared with no balcony: 


=f) 20 —(—5) 2 

— 2 : — — OF 
wy ff Los fo 
+ +— + 


h, k, h, 8 1 20 


The difference for one balcony is A, (77.2 — 52.6) = 0.8 x 24.6 =19.7W 
For 350 apartments, the difference is 6897 W. 


For the steel supported balcony where A, = 0.01 mand P=0.6m 
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As before, however, in this case Bi << J because k, >> k, 


1/2 1/2 
mW = ae w= zu x2211 
kA 40x 0.1 


mW >> 2, so we can use the infinite fin approximation as before 


ane sas 
(ee cal =182 Wim? 


ttf AY O34 0.01 ) 
he kee 8 40 | 20x6x40 


L 


O, = 4.q, =9.01x182 =1.82 W/beam 


For 350 apartments, O, =1915 W 


Example 2.10 
In free convection, the heat transfer coefficient varies with the surface to fluid temperature difference 
(7 only ). Using the low Biot number approximation and assuming this variation to be of the form 


h= G(r. te )" Where G and 7 are constants, show that the variation of the dimensionless 


temperature ratio with time will be given by 


O-" =1+(nh,,,A)t 


init 


Where 
= (7, -T a) 7 Area 
Tina — Ty Mass x Specific Heat Capacity 
and h,,,, = the heat transfer coefficient at t= 0. Use this expression to determine the time taken for an 


aluminium motorcycle fin (9 = 2750 kg/m*, C = 870 J/kgK ) of effective area 0.04 m’ and 


thickness 2mm to cool from 120°C to 40°C in surrounding air at 20°C when the initial external heat 


transfer coefficient due to laminar free convection is 16 W/m* K. Compare this with the time estimated 


~") which assumes a constant value of heat transfer coefficient. 


from the equation (@ = e 
Solution 


Low Biot number approximation for free convection for B, << 1 


Heat transfer by convection = rate of change of internal energy 
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d(T, —T,) 
hA(T, ~T,) =-mC——— (1) 


We know that h = G(T, —T,)" 


Where G is a constant. 


(Note that this relation arises from the usual Nusselt/Grashof relationship in free convection; for 


example: Nu = 0. \(Gr Pr)” in turbulent flow or Nu = 0.54(Gr Pr)'’* for laminar flow) 
Equation | then becomes: 


A dt 


GP a7, (= t= 


(ot d(T, -T,) 


2, mC 7 t-0 (7, 7 T, ye 
GnA a = 
met, -7,)" -@,-2, Q) 


At t= 0, (7, -7,)=(Z,,-T,) 


If we divide equation 2 by (r ery i 


And use the definition 0 = (7, -T,) 
We obtain 
GnAt GnAt 
= Qo” l= T n 
mC(L, —T, y mC ( - ') 


Since G(r. ae ri f ) =h,, the heat transfer coefficient at time t = 0, then 


h,At 
= + 
mC 


go” 1 
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Conduction 


Or 6" =nhAt+1 


h=16W/m’?>K 
eg ay 


For aluminium p = 2750 kg/m’, 


For laminar free convection, n= % 


m= pAX =2750x0.04x 0.002 = 0. 
is Fay 
mC 0.22x870 


0” =nh,At+1 which gives 
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nh,A 
When T=40°C 9 = 20720 _ 
120 —20 
Then 
“1/4 
_ (0.21) seme 


(1/4)x16x 2.1x10~ 


For the equation 0 =e" 


which assumes that the heat transfer coefficient is independent of surface-to-fluid temperature 
difference. 


Ind In 0.2 


= = = 4795 
—-hA —-16x2.1x10% 


t 


Percentage error = “_ x 100 =19% 


Example 2.11 


A 1 mm diameter spherical thermocouple bead (C = 400 J/kg K, p = 7800 kg m?) is required to 
respond to 99.5% change of the surrounding air p = 1.2kgm?, w=1.8x10%kg/ ms , k= 
0.0262 W/m K and Pr = 0.77) temperature in 10 ms. What is the minimum air speed at which this 
will occur? 


_—- d=3mm 


Download free ebooks at bookboon.com 


32 


Heat Transfer: Exercises Conduction 


Solution 
Spherical bead: area = nd? 
volume = md?/6 


Assume this behaves as a lumped mass, then 


Step change 
In air 
temperature 


Time, t 
Ty 3 
= 0.995 
Too — T; 
(given) 


For lumped mass on cooling from temperature 7; 


ES 1 (—A t) = 0.995 
T,-7,  °%P = 0. 
8 t = 0.01 

> ae 7 =U. S 
At = 0.005 

7208 


Which gives the required value of heat transfer coefficient 


hA 


——- = 05 
pvc 


So 


Download free ebooks at bookboon.com 


33 


Heat Transfer: Exercises Conduction 


Aes 0.5 x 1073 x 400 x 7800 


- = 260 W/m? K 


vy, = MD 260X107? _ 
es Oe WT 


For a sphere 


Nup = 2 + {0.4Rep!? + 0.06Re;/*} Pr°4 


From which with Pr = 0.707 
f =0.4Re,/” + 0.06Re2/? -9.4 =0 
f = 0.2Re,/* + 0.04Re, ”? 


Using Newton iteration 


x tt) = yr f(x) 


f(x) 


Starting with Rep = 300 


0.4V¥300 + 0.06(300 2/3 _g 4 0.222 
Ree? sng = [0.4V300 + 0.06(300)/" — 9.4] “—39h0 = 
D 0.2 0.04 | 0.01782 
¥300 3001/3 


Which is close enough to 300 


From which 
Re 
eS eeiae 4.5 m/s 
Dp 
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3. Convection 


Example 3.1 


Calculate the Prandtl number (Pr = wC,/k) for the following 


a) 
b) 
c) 


d) 


e) 


g) 


Water at 20°C: w= 1.002 x 10° kg/ms, C, = 4.183 kJ/kg K and k = 0.603 W/m K 

Water at 90°C: p= 965 kg/m’, v = 3.22 x 10°’ m’/s, C, = 4208 J/kg K and k = 0.676 W/m K 
Air at 20°C and 1 bar: R = 287 J/kg K, v = 1.563 x 10° m’/s, C, = 1005 J/kg K and 

k = 0.02624 W/m K 

_ 1.46x10°T*” 
~ (110+7) 


C, = 0.917 + 2.58x 10*7T —3.98x10°T* kJ/kg K (Where T is the absolute temperature in 


K) and k=0.03186 W/m K. 

Mercury at 20°C: w= 1520 x 10° kg/ms, C, = 0.139 kJ/kg K and k = 0.0081 kW/m K 
Liquid Sodium at 400 K: = 420 x 10° kg/m s, C, = 1369 J/kg K and k= 86 W/mK 
Engine Oil at 60°C: w= 8.36 x 10” kg/m s, C, = 2035 J/kg K and k= 0.141 Wim K 


Air at 100°C: kg/m s 


lly thanks to our 
nd automatic 
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Solution 
Cc 3 
a) Prat 7 _ 1.002 x10 x 4183 _ ¢ gs 
k 0.603 
C vC 7 
b) peal -aF p _ 9653.22 x10 x 4208 _ 1 93 
k k 0.676 
C 
c) Pr= aie 
k 
P _ 100000 


=1.19 kg/m? 


PRT 287x293 


-5 
_ 1.19x1.563x107 x1005 0.712 


Pr 
0.02624 


~6 4793/2 6 3/2 
_146x10° 7? _1.46%10° x373"? _ 9 19105 ke/ms 
(110+7) 110 +373 


d) 


C, = 0.917+2.58x107T —3.98x10°T* = 0.917 + 2.58x10™ x 373 —3.98x10™ x373" 
= 1007.7 J/kg K 


-5 
a” 2.18x10~° x 1007.7 ~ 0.689 


0.03186 


Pr 


C 6 
HU, _ 1520x110 x19? 0.0261 
k 0.0081 x 10 


e) Pr= 
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MC, 420x10° x1369 


f) Pr= = 0.0067 
k 86 
Cc 8.36x107 x 2035 
epee ee ee) 
k 0.141 
Comments: 
e Large temperature dependence for water as in a) and b); 
e small temperature dependence for air as in c) and d); 
e use of Sutherland’s law for viscosity as in part d); 
e difference between liquid metal and oil as in e), f) and g); 
e units of kW/m K for thermal conductivity; 
e use of temperature dependence of c, as in part a). 
Example 3.2 


Calculate the appropriate Reynolds numbers and state if the flow is laminar or turbulent for 


the following: 


a) 


b) 


c) 


d) 


) 


A 10 m (water line length) long yacht sailing at 13 km/h in seawater p = 1000 kg/m’ and 

H=13x 10° kg/ms, 

A compressor disc of radius 0.3 m rotating at 15000 rev/min in air at 5 bar and 400°C and 
_ 1.46x10°T*” 
~ (110+7) 

0.05 kg/s of carbon dioxide gas at 400 K flowing in a 20 mm diameter pipe. For the viscosity 

take Ww = 1.56x10°7" | om S 

(233+T) 


The roof of a coach 6 m long, travelling at 100 km/hr in air (9 = 1.2 kg/m’ and w= 1.8 x 10° 


kg/m s 


kg/m s) 

The flow of exhaust gas (p = 1.1 bar, T= 500°C, R = 287 J/kg K and = 3.56 x 10° kg/m s) over 
a valve guide of diameter 10 mm in a 1.6 litre, four cylinder four stroke engine running at 

3000 rev/min (assume 100% volumetric efficiency an inlet density of 1.2 kg/m? and an exhaust 
port diameter of 25 mm) 
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Solution 
;_ 13x10° sib 
i 
). Roe Cee =2.78x107 (turbulent) 
L 3x 


b) 7 =400+273=673 K 


a6 3/2 
_ 146x107 x673") _ 3 56,1055 kg/m s 
(110+ 673) 


\ 
15000 
= ———_ x 

60 


@) 22 =1571 rad/s 


u=Qr=1571x0.3 = 471.3 m/s 


_ P _ 100000 
PRT 287x673 


= 2.59 kg/m?* 
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Characteristic length is r not D 


Rea ee (turbulent) 
uu 3.26% 107 
2 
c) m= puA=pux 4 
An 
paD* 
m 
Re = PUP _ px4mD _ 4m sd 
MH pmD?y mDy 


_ 1.56x10% x 400°” 
(233 + 400) 


=1.97x10° kg/ms 


Re = pac =1.6x10° (turbulent) 


mx 0.02 x1.97x10> 


— 100x10° 
3600 


d) =27.8m/ 5 


— pul _1.2x27.8x6 


aa =11.1«10’ (turbulent) 
u 8x 


Re 


e) Let m be the mass flow through the exhaust port 


m= inlet density X volume of air used in each cylinder per d 1 p D 
second 


1.6x10° 3600 1 
- x _ 


m=1,.2 = 0.012 kg/s 
60 2 
_ 4m 
aD’ p 
Re =e 
ya 


Download free ebooks at bookboon.com 


39 


Heat Transfer: Exercises 


Convection 


— 4x 0.01 au = 6869 (laminar) 
ax3.56x107 x 0.025 


Comments: 


e Note the use of D to obtain the mass flow rate from continuity, but the use of d for the 
characteristic length 


e Note the different criteria for transition from laminar flow (e.g. for a pipe Re ~ 2300 plate 
Re x 3x10°) 


Example 3.3 


Calculate the appropriate Grashof numbers and state if the flow is laminar or turbulent for the following: 


a) 


b) 


c) 


d) 


A central heating radiator, 0.6 m high with a surface temperature of 75°C in a room at 18°C (o= 
1.2 kg/m? , Pr = 0.72 and w= 1.8 x 10° kg/m s)] 

A horizontal oil sump, with a surface temperature of 40°C, 0.4 m long and 0.2 m wide containing 
oil at 75°C (9 = 854 kg/m’ , Pr = 546, B=0.7x 10° K' and w=3.56 x 10° kg/m s) 

The external surface of a heating coil, 30 mm diameter, having a surface temperature of 80°C in 
water at 20°C (p= 1000 kg/m’, Pr = 6.95, B= 0.227 x 10 °K! and w= 1.00 x 10°kg/m s) 

Air at 20°C (p= 1.2 kg/m’ , Pr = 0.72 and z= 1.8 x 10° kg/m s) adjacent to a 60 mm diameter 
vertical, light bulb with a surface temperature of 90°C 


Solution 


a) 


b) 


29 BATL 
Geena 


a 
AT =75-18=57K 


pale ae 
T (18+273) 291 


_ 1.2? x9.81x 57x 0.6° 


Gr —— =1.84x10" 
291x(1.8x10°) 
Gr Pr =1.84x10" x 0.72 =1.3x 10° (mostly laminar) 
a oma _ 0.4x 0.2 ~ 0.0667 m Oil 
Perimeter 2.x (0.4 + 0.2) g °C 
T, = 40°C 
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AT =75—-40 =35K 


Ge pgeBATL 854? x9.81x0.7x10~ x35 x 0.0667" 
2 : = 


fe (356x107) 
Gr Pr = 4.1x10* x 546 = 2.24x10’ 


=4.1x10' 


Heated surface facing downward results in stable laminar flow for all Gr Pr 


c) 
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d) 


AT = 80—20=60K 


—_ p gBATL 1000? x9.81x0.227x 10% x 60x 0.03° 


2 x10) = 3.6x10° 


Gr Pr = 3.6 x10° x 6.95 = 25x10° (laminar) 


ya Area _mD°? _D a 
Perimeter 47D 4 2 


AT =90-20=70K 


1 ae ae ) \ 
P= == 04273) 293 — 


Re pg BATL _1.27x98x70x0.015° 


2 


Ml 293 x(1.8x10%) 


=3.5x104 


Gr Pr =3.5x10* x 0.72 = 2.5x10* (laminar) 


Comments: 


e Note evaluation of # fora gas is given by 6 =1/T 


e Fora horizontal surface L = A/p 


Example 3.4 


Calculate the Nusselt numbers for the following: 


a) 


b) 


oy) 


d) 


A flow of gas (Pr = 0.71, = 4.63 x 10° kg/m s and C, = 1175 J/kg K) over a turbine blade of 
chord length 20 mm, where the average heat transfer coefficient is 1000 W/m’ K. 

A horizontal electronics component with a surface temperature of 35°C, 5 mm wide and 10 mm 
long, dissipating 0.1 W by free convection from one side into air where the temperature is 20°C 
and k= 0.026 W/m K. 

A 1 kW central heating radiator 1.5 m long and 0.6 m high with a surface temperature of 80°C 
dissipating heat by radiation and convection into a room at 20°C («4 = 0.026 W/m K assume black 
body radiation and o= 56.7 x 10’ W/m K*) 

Air at 4°C (4 = 0.024 W/m K) adjacent to a wall 3 m high and 0.15 m thick made of brick with k= 
0.3 W/m K, the inside temperature of the wall is 18°C, the outside wall temperature 12°C 
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Solution 
C 
a) Pr= se 
k 
C -5 
peor p _ 4.63x10 x75 _ 9 9766 W/m K 
Pr 0.71 
ae AL _ 1000 x 0.02 6] 
k 0.0766 
b) Nu _AL_@t 
k AT k 
jo2= os = 2000 W/m? 
A 0.01x 0.005 


AT = 35-20 =15°C 
_ Area 50 5 


- mm = 0.001667m 
Perimeter 30 3 


AL _ 2000x0.001667 _ 


Nu 8.5 
k 15x 0.026 
q, L 
c) Nu=t— 
AT k 


In this case, q must be the convective heat flux — radiative heat flux 


T, =80+273=353 K 

T,, = 20+273=293 K 

O, = 0A(T’ -T!)= 56.7x10~ x (1.5 0.63534 - 293')= 416 W 

AT =80-20=60K se me! 
' 


O. =O-O, =1000—416 =584W 


Tr, - 80°C 


Download free ebooks at bookboon.com 


43 


Please click the advert 


Heat Transfer: Exercises 


Convection 


q. Mei a 8 = 649 W/m? 
A 1.5x0.6 
ieee q. L_ 649 0.6 — 249 
ATk 60 0.026 
d) AT=12-4=8K 
g= abe =) eae 
W 
(assuming 1-D conduction) 
_ O38 =12) _ spay 
0.18 
Miata 2188 
ATk 8 0.024 
Comments: 


k,=0.3W/mK 


I, 


|-—___4 


W=0.15m 


e Nuis based on convective heat flux; sometimes the contribution of radiation can be significant 


and must be allowed for. 


e The value of k is the definition of Nu is the fluid (not solid surface property). 


e Use of appropriate boundary layer growth that characterises length scale. 


What do you see? 


A skateboarder? 
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Example 3.5 


In forced convection for flow over a flat plate, the local Nusselt number can be represented by the 


general expression Nu, = C, Re”. In free convection from a vertical surface the local Nusselt number 


is represented by Nu, = C,Gr,” , where C), C2, n and m are constants 


a) Show that the local heat transfer coefficient is independent of the surface to air temperature 
difference in forced convection, whereas in free convection, 4, depends upon (7, — T..)” 


b) In turbulent free convection, it is generally recognised that m = 1/3. Show that the local heat 
transfer coefficient does not vary with coordinate x. 


Solution 
a) Nu, = ua 
k 
Re, = — 
L 


Hence Ko (22) 


This shows that the heat transfer coefficient for forced does not depend on temperature difference. 


For free convection Nu, =C,Gr" 


2 3 
Gr 2 SBATX 


x uw 
2 3\" 
Hence h= Ao (2 anaes (1) 
x Hu 


So for free convection, heat transfer coefficient depends on AT” 


b) From (1), with m = 1/3 for turbulent free convection: 
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2 


2 3 1/3 
nko? g BAT x 
x Ul 


2 


k 4 BAT 1/3 
h= o[ #8 x 
vi 


2 BAT 1/3 
p=ace| PEAT ae 
vi 


Hence the convective heat transfer coefficient does not depend on x 
Example 3.6 


An electrically heated thin foil of length L = 25 mm and width W = 8 mm is to be used as a wind 


speed metre. Wind with a temperature 7), and velocity U,, blows parallel to the longest side. The foil 

is internally heated by an electric heater dissipating O (Watts) from both sides and is to be operated in 
air with T,, = 20°C, C, =1.005kJ /kg K, v =1.522x10°m’/s p=1.19kg/m? and Pr =0.72 
. The surface temperature, 7°, of the foil is to be measured at the trailing edge — but can be assumed to 


be constant. Estimate the wind speed when 7., = 32°C and O =0.5W. 


. 
u,, 

Solution ——> 

. & 
Firstly, we need to — ce) 
estimate if the flow iS 
laminar or re 
turbulent. 25 mm 


Assuming a critical (transition) Reynolds number of Re = 3 x 10° the velocity required would be: 


_3x10°3x10%v _ 3x10°x1.522x10° 
me pe - i 25x107 


u = 304 m/s 


Wind speed is very unlikely to reach this critical velocity, so the flow can be assumed to be laminar. 


Nu, =0.331Re\? Pr’ 


Yay 


Nu.,, = 0.662 Re)” Pr!’? = (rT. -T.k 
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Jab 


Rel’? _ 
"  (T.-T,,)\k x 0.662 Pr’? 


~__ 05/2 _ 1959 W/m? 


” 0.025 x 0.008 


Rell? = 1250 x 0.025 173.5 
(32 —20)x 0.0253 x 0.662 x (0.72) 


Re, =3x10* 


= Ren¥ _ 3x10! x1.522x10% 


= = ne =18.3 m/s 
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Example 3.7 


The side of a building of height H = 7 m and length W = 30 m is made entirely of glass. Estimate the 
heat loss through this glass (ignore the thermal resistance of the glass) when the temperature of the air 
inside the building is 20°C, the outside air temperature is -15°C and a wind of 15 m/s blows parallel to 
the side of the building. Select the appropriate correlations from those listed below of /ocal Nusselt 
numbers to estimate the average heat transfer coefficients. For air take: p= 1.2 kg/m’, = 1.8 x 10° 
kg/ms, C, = 1 kJ/kg K and Pr = 0.7. 

e Free convection in air, laminar (Gr, < 10”): Nu, = 0.3 Gr," 
e Free convection in air, turbulent (Gr, > 10°): Nu, = 0.09 Gr," 
e Forced convection, laminar (Re, < 10°): Nux = 0.33 Re,”? Pr’? 
e Forced convection, turbulent (Re, > 10°): Nu, = 0.029 Re,”* Pr!” 


inside 
T, 
‘ outside 
s 
h, 
u, =15 m/s 
Solution 
Cc Cc 1.8x10° x1 
pe gives: a= Bias DOES 000 | 0.026 W/mK 
k Pr 0.7 


First we need to determine if these flows are laminar or turbulent. 
For the inside (Free convection): 


get-_! 1 og 
—T (204273) 293 
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pie pgBATL 1.2°x981xATx7 


2 


a (18x10) x293 


Gr =5.1x10'° AT 
(Flow will be turbulent over most of the surface for all reasonable values of AT ) 
For the outside (Forced convection) 


_ pub —1.2x15x30 
1.8x10° 


Re, =3x10’ 


(Flow will be turbulent for most of the surface apart from the first 0.3 m) 


Hence we use the following correlations: 


On the inside surface: Nu, =0.09Gr"" 


On the outside surface: Nu, = 0.029 Re?* Pr'’* 


For the inside: 


2 _ 3 1/3 
ee -4X-o02 g BT. Ze) 


2; 
a 
3\13 
x 
h =constant x () = constant 
x 


Hence heat transfer coefficient is not a function of x 


i =f; (1) 
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x=L x=L h 
ee [ade = . [xPdx = = (2) 
Ly) OLS, 0.8 


Write a heat balance: 


Assuming one-dimensional heat flow and neglecting the thermal resistance of the glass 


q=h,(T, -T,) 
q=h,(T, -T,) 
h(T,-T,)=h,(T, -T,) (3) 
From equation 1 
HE ate pial, T.)H\ 
k wOXE 
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2 — 
h, =0.09| 2 eats 1) | ,.0.026 
(18x10) x293 
h, =1.24(7, -7,)"” (4) 


From equation 2: 


0.8 
h,W _ 0.029( puW pr!3 
k 0.8 a 


0.8 
jy = 02026 , 0.029 (1.21530) 9 ais 
° 300.8 \ 1.810 
h, = 26.7W/m? K (5) 


From (3) with (4) and (5) 

1.24(7, -T,)* =26.7(T, -T,) 

1.24(20-T,)"* =26.7(T, +15) 

T, = 0.0464 (20-7, )** -15 (6) 
To solve this equation for 7’, an iterative approach can be used 

First guess: T, =—10°C 

Substitute this on the right hand side of equation 6: 

T, = 0.0464 (20-—(-10))"* -15 =-10.7°C 

For the second iteration we use the result of the first iteration: 

T, = 0.0464 (20-—(-10.7))"* —15 =-10.6°C 


The difference between the last two iterations is 0.1°C , so we can consider this converged. 


T, x -10.6°C 
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From which: 
q=h,(T, -T,)=26.7(-10.64+15)=117 W/m? 
O = gA =117x30x7 = 24600 W = 24.6 kW 


Example 3.8 


The figure below shows part of a heat exchanger tube. Hot water flows through the 20 mm diameter 
tube and is cooled by fins which are positioned with their longest side vertical. The fins exchange heat 
by convection to the surrounds that are at 27°C. 


Estimate the convective heat loss per fin for the following conditions. You may ignore the contribution 
and effect of the cut-out for the tube on the flow and heat transfer. 


a) natural convection, with an average fin surface temperature of 47°C; 
b) forced convection with an air flow of 15 m/s blowing parallel to the shortest side of the fin and 
with an average fin surface temperature of 37°C. 
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The following correlations may be used without proof, although you must give reasons in support of 
your choice in the answer. 


Nu, = 0.3 Re’? Pr’? Re, <3 x 10° 
Nu, = 0.02 Re,”* Pr’? Re, >3 x 10° 
Nu, = 0.5 Gry! Pr! Gr, < 10° 
Nu, = 0.1 Gr,!? Pr? Gr, > 10° 


For air at these conditions, take: Pr = 0.7, k = 0.02 W/mK, p= 1.8 x 10° kg /ms and p= 1.0 kg/m? 


Solution 


100 mm 


20 mm dia 


On the outside of the water tube, natural convections means that we need to evaluate Gr number to see 
if flow is laminar ot turbulent 


2 3 
ATL 
Gk & Pi 


Wi 
AT =47-27=20K 


1 La 


P= TT 973 300 


2 3 
Gra! x 9.81x 20x 0.1 ~2x108 (Laminar) 


(1.8x10°} «300 


(L here is height because it is in the direction of the free convection boundary layer) 
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So we use: 


Nu, =0.5(Gr, Pr)" 


i L. 
i= +f nde = constant [a ax 
L 0 0 


2 


a 


Nu,, = —(2x10° x 0.7)" =23 


Nu,,k _ 23x 0.02 
a | 0.1 


=46W/m’> K 


Day = h,, AT 
O= ¢,A=h,, ATA =4.6x20x0.1x0.05x2 (Last factor of 2 is for both sides) 


Q=0.92 W 


For forced convection, we need to evaluate Re to see if flow is laminar or turbulent 


S 
ey 


100 mm 


20 mm dia 
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— pul _1x15x0.05 


i8xi07 = 4.17x10* (Laminar) 
u 8x 


Re 


(L here is the width because flow is along that direction) 


Nu, =0.3Rel? Pr!” 


Nu,, = 0.6Rel? Pr!’ =0.6x (4.1710) |” x 0.7! = 109 


p= Nak _109%002 _ 45 5 yy jm? 
L 0.05 

OQ =4,,A=h,,ATA = 43.5x10x0.10.05x2 (AT =10°C) 

0=435W 

Example 3.9 


Consider the case of a laminar boundary layer in external forced convection undergoing transition to a 
turbulent boundary layer. For a constant fluid to wall temperature difference, the /ocal Nusselt 


numbers are given by: 


Nu, = 0.3 Re,’ Pr’? (Rex < 10°) 
Nu, = 0.04 Re,”* Pr'? (Re, > 10°) 


Show that for a plate of length, L, the average Nusselt number is: 


Nuay = (0.05 Re,°* - 310) Pr’? 
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Solution 


Where for a constant surface-to-fluid temperature: 


x L 
hiz= - J Pasina + Laer dx 
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Since for laminar flow (Re, < 10°): 


Nu, =0.3Rel? Pr!” 


h 


lam 


1/2 
=0.3x ie) Pe ag SC a? 
u 


Where C,, does not depend on x 


lam 
Similarly: 


~ -0.2 
a Cus x 


0.8 
Cs = 0.04~x i Pe Pr’? 


HM 
Hence 
l xX, L 
ae -1/2 ~02 
hy, = L J Cio dx a ! Ciyh dx 
1 1/2] os |4 
a art Cram + Cun uae 
L) ae. 0.8 |. 
hk 
Nu, = 
L 
Cc Cc 
Nu, — ~lam oe eg Ea turb [2° — x8] 
k 0.8k 


1/2 0.8 0.8 
Mu, = a{ Pt xt? Pr'/34.0.05 wa [ea Pr 
Le i i 
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=10° (The transition Reynolds number) 


So 

Nu,, = Pr’? 0.6 x (10°)? +.0.05Re?*— 0.05 x (10° | 

Nu,, =(0.05Re°*-310)Pr! 

Example 3.10 

A printed circuit board dissipates 100 W from one side over an area 0.3m by 0.2m. A fan is used to 
cool this board with a flow speed of 12 m/s parallel to the longest dimension of the board. Using the 
average Nusselt number relationship given in Example 3.9 to this question, calculate the surface 
temperature of the board for an air temperature of 30 °C. 


Take an ambient pressure of | bar, R = 287 J/kg K, 


Cp = 1kJ/kg K, k= 0.03 W/m K and p=2 x 10° kg/ms 


0.2m 


q=100W 


u, =12m/s 


Ir =30C 
Solution 
da = 2 a ane = 1666.7 W/m? 
A 0.2x0.3 
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— uC, — 2x10 x10° 
k 0.03 


Pr = 0.667 


P 10° 


= = =1.15kg/m 
RT 287x303 


ye) 


dels x12x0.3 


Re, = =—=2.07x 10° 
2x10 


Using the formula for Nusselt Number obtained in Example 3.9: 


Nu,, = (0.05 Reo*—310)Pr'”? 


Nut, = (0.05 x 2.07x10-5}"* - 310} (0.667)'" =511 


hk dak 
Nu,, = av =_ Day 
L ATk 
pep GaP SNOCOTROS 5 poe 
Nu,.k 511x0.03 
T. =T,+AT 


T, = 30+ 32.6 = 62.6°C 
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4. Radiation 


Example 4.1 


In a boiler, heat is radiated from the burning fuel bed to the side walls and the boiler tubes at the top. 
The temperatures of the fuel and the tubes are 7, and T> respectively and their areas are A; and Ap. 


a) Assuming that the side walls (denoted by the subscript 3) are perfectly insulated show that the 
temperature of the side walls is given by: 


T, = [Af + A,F,,T,' ) 
ie 
A,F,, + AF; 


where F)3 and F; are the appropriate view factors. 


b) Show that the total radiative heat transfer to the tubes, Q,, is given by: 


A, F,,4,F. 
Q, = [ar Athol -T;) 


A,F, + AF 


6007 'A'G Swewsds Wy] 222] @ 


iB \ 
3 It’s only an 


opportunity if 
you act on it 


IKEA.SE/STUDENT 


eN 10 7 
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c) Calculate the radiative heat transfer to the tubes if 7, = 1700°C, T, = 300°C, A; = A2 = 12m? and 
the view factors are each 0.5? 


tubes (2) 


Adiabatic wall (3) 


Fuel 1) 


Solution 

a) O,=0,,+03, (1) 
Since the walls are adiabatic 
Os. =O. (2) 
From (2) 
o A, Fy(T4 -T3)=0 4, F(T -T/) 


_A Fig Ty + Ay Fy Ty 


i 
A, Fy + A, Fi, 


1/4 
since A, F, =A, F 


T. [AREA 
, = 
A, F,, + A, Fi 


b) From (1) 


QO, =0A4, Pole -Ti)+oA, Fi \t -T;) 


QO, =0A4, Fe -Tj)+oA, PAT -T;) 
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QO, =o A, ele TF Ga F,, 


AFT, tA r | 
2 
A, Fy, +A, F,, 


QO, =o A, F(T -T;)+0A, Fi; 


A, Fy + A, Fi; 


(4 Fi; a + A, Fy, ie — A, ony i — A, 1) 
Q, =o A, F(T; -T;)+oA, P| 


A, F,, ies =A om) 
A, Fy +A, F, 


i A F 
0, =0 A F(T) -T;)+o A, F(T; ri] — 


A, Fy + A, Fi 
: A, F,,A, F, 
0, =0 4, F,(T' -T})+o(r) ii 2 £93441 £13 
A, Fy + A, Fi 
. A, F,,A, F 
Q, =o(T' -T; AF, +2 
A, F,, + A Fi 


_A FT! +4; Fy T _A FT) +4, FT, 


c ee 
° A, Poy + A, #); Ay PPA Ps 


_ 12x 0.5x1973* +12x0.5x 573" 
12x 0.5+12x0.5 


T; = 1662 K 


QO, = 56.710” (19734 -s73*[6+°8) = 7.68x10° W 
+ 


Example 4.2 


Two adjacent compressor discs (Surfaces | and 2) each of 0.4 m diameter are bounded at the periphery 
by a 0.1 wide shroud (Surface 3). 


a) Given that F), = 0.6, calculate all the other view factors for this configuration. 
b) The emissivity and temperature of Surfaces | and 2 are ¢; = 0.4, T; = 800 K, & = 0.3, T2 = 700K 
and Surface 3 can be treated as radiatively black with a temperature of T; = 900 K. Apply a grey 


body radiation analysis to Surface | and to Surface 2 and show that: 


2.5 J, -0.9 Jn = 45545 W/m? 
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and 
3.333 Jo—1.4J3,=48334 W/m’. 


The following equation may be used without proof: 


(a Aa: 
— =S F J,-J, 
ee ae 


E. 


L 


c) Determine the radiative heat flux to Surface 2 


é,=0.4 
T, =800K 


Solution 


a) r=r,=r=02m 


a=0.lm 
rm _02_, 
a 0.1 
Eo agg 
r, 0.2 


Ff, =90.6 (Although this is given in the question, it can be obtained from appropriate tables 


with the above parameters) 
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Ff, =9 (As surface | is flat, it cannot see itself) 
Ff, =1-0.6=0.4 (From the relation ye = | in an enclosure) 


Fy, =0.6 (Symmetry) 


Fy =0 
F,, =0.4 
2 
jee ie hha 


A, 8 2xax0.2x01 


F,, =0.4 (Symmetry) 


F,, =1-0.4-0.4=0.2 
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Radiation 


b) 


é; 
(E,, -J,)=¢[F. VU, -42)+ Fa - Js) 
Ey, =F {1+ OF + O,Fia}- OF ode - OF Js 
Bg =67, 


J,=o0T; (Radiatively black surface) 


oT’ =2.5J,-0.9J, -0.60T; 


56.7x10° x 800* = 2.5x J, -0.9x J, -0.6 56.7 x 107 


2.5J, -0.9J, =45545 W/m? 

Applying to surface2 (i=2) 

E,o =I 14+ OF + OF} FJ) — bo FJ 
La=oT, 


l=—e, _ 1-03 
Ey 0.3 


= 2.333 


p, = 


oT, =3.333 J, -1.4J, -—0.93330T; 


3.333 J, -1.4J, = 48334 W/m? 


> «9004 


(1) 


(2) 
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c) From (2): 


_ 3.333. J, — 48334 
1.4 


J 


Substituting in (1) 


: 3.333 J, — 48334 
1.4 


2.5 0.97, = 45545 W/m? 


J, = 26099 W /m? 
The net radiative flux to surface 2 is given by 


_ Ex ,-J, _ 56.710” x 7004 — 26099 
a ae 1-03 
£5 0.3 


= —5.351x10°? W/m? 


The minus sign indicates a net influx of radiative transfer as would be expected from 
consideration of surface temperatures. 
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Example 4.3 


The figure below shows a simplified representation of gas flame inside a burner unit. The gas flame is 
modelled as a cylinder of radius r, = 10 mm (Surface 1). The burner comprises Surface 2 (a cylinder of 
radius r, = 40 mm and height h = 40 mm), concentric with Surface 1 and a concentric base (Surface 3), 
of radius r; = 40 mm. The end of the cylinder, Surface 4, opposite to the base is open to the 
surrounding environment. 


a) Given that F,; = 0.143 and F,. = 0.445 use the dimensions indicated on the diagram to calculate 
all the other relevant view factors. 


b) The flame, base and surroundings can be represented as black bodies at constant temperatures 77, 
T; and T,, respectively. The emissivity of the inside of Surface 2 is e€, = 0.5. Apply a grey body 
radiation analysis to Surface 2 and show that the radiosity is given by: 


fam o(T; +F,,T;' +F,,T, +F,,T,') 
; [ey eg PP a 
The following equation may be used without proof: 


c) The temperatures T, and T3 are found to be: T; = 1800K and T; = 1200K, and the surrounds are at 
500 K. Estimate the temperature T>, using a radiative heat balance on the outer surface of Surface 
2, where the emissivity is €9 = 0.8 


r, = 40mm 


¥,=10mm 


Surface 1 


ae 


h=40mm 


_L 


Surface 4 


Surface 2 


Surface 3 


Solution 


a) A, =2anh 
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A, =2a19,h 


A, =A, =a(? -7?) 


Borg tliat hy =] 
but 
AF, > AF, 


F,=“2F, =F, =) x0.14338 = 0.57352 
A r 7 10 


1 1 


Thus 


F.=F _ 1-0.57352 _ 9 54394 
13 14 5 


Fy, + Fy + Py t+ Fy =1 


F, =F, = 1 Fy Fy _1-0.14338-0.44515 _) aneny 
5) 2 


F,, =0 
A Fi; — A;F3, 
Bop eg oo 011 


A? ae)? 0.04 -0.07 


A,F,, = AF, 
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Radiation 


Fy, =1-—0.11373 — 0.43891 = 0.44736 
Similarly (using symmetry) 


F,, = Fy, = 0.11373 


Fy = Fy = 0.43891 


Fy = Fy, = 0.44736 


F,, =0 
i ae 
é 


69 


_ 2x0.04x 0.04 
A, a(re-r?) * 0.047 -0.07 


x 0.20574 = 0.43891 
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Radiation 


c) 


For surface 2,1= 2, j=1,3,4 


E,,-J 
a = Fala — Si) 4 Fas U2 - Is) + Faso - I) 
2 
&4 
1-0. 
é€, =0.5, La 
0.5 
J, =E,,, J3=E,3, J, =F,, (1,3, 4 are black) 


Fy =d5 Fis = Sa OP SE ei Ar i) 
JF, + Fy t+, tl=oT, +oT{F,, +07, F,, +07, F,, 


o(T! 47°F, 47 F,, +11 F,) 
Po tiie Phy +l 


J, = 


_ 56.7x10° (T4 +18004 x 0.57352 +1200 x 0.20574 +500* x 0.20574) 


a 
0.57352 + 20574 + 0.20574+1 


J, =36.47x10°T,' + 70913 
On the outside of surface 2: 
-q,= oe, (Te 27) 

Also 


E,,-J 
2 = oT; -36.47x10°T;' + 70913 


qy = l-€é, 


é> 
20.23x10°T;' + 70913 = -56.7 x10” x 0.8(7, - 500°) 


T, =1029K 
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Example 4.4 


The figure below shows a schematic diagram, at a particular instant of the engine cycle, of a cylinder 
head (Surface 1), piston crown (Surface 2) and cylinder liner (Surface 3). 


a) Using the dimensions indicated on the diagram, and given that F). = 0.6, calculate all the other 
relevant view factors. 


b) The cylinder head can be represented as a black body at a temperature T; = 1700 K and the 
emissivity of the piston crown is€é, = 0.75. Apply a grey body radiation analysis to the piston 


crown (Surface 2) and show that the radiosity is given by: 
Jp = 42.5 x 10° T:* + 71035 + 0.1 J; 
The following equation may be used without proof: 


Eyad: 
Cae 
i i J= 


c) Similar analysis applied to the cylinder liner gives: 
J; = 107210 + 0.222 Jp 


If the surface temperature of the piston crown is, 7, = 600 K, calculate the radiative heat flux into 
the piston crown. 


d) Briefly explain how this analysis could be extended to make it more realistic 


Surface 1 — 100 mm — 
Cylinder Head 


(1700 K) 
Surface 3 25mm 
Cylinder liner 


¥ 


Surface 2 
Cylinder piston 
crown (600 K) 


Solution 


a) A =A, =ar° =2x50° =25002 mm’ 
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Radiation 


A, = DL = 7 x100x25 = 25002 mm? 


F,, =0 (Flat surface) 


fF, =0.6 (Given) 


F, =1.0-F,, =1 


By Symmetry: 


F,, =F, = 0.6 
F,, =F, =0A4 
F,, =0 

A, 


Fy =F =04 Since 4, = 


3 


.0—-0.6 =0.4 
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Radiation 


b) 


F,, =0.4 (By symmetry) 
F,, =1.0-—F,, —F,, =1.0-0.4-0.4=0.2 
For surface 2, 1=2 


E,o -J, 


= = Fi (In Ji) alas) 
a2 


a) 


J, =oT,' (Black body) 


ey ee 
0.75 3 
Lg=ol, 
el. ag. 
alla et er: —J;) 


1 
ois 3 fue! +F,,J;) 


J, = 
1+ 3% + F,;) 


J,= 


56.710 xT; + (0.636. 10° x1700* +0.4J; ) 


1+ 5(0.6+04) 


J, =42.5x10°T; + 71035 +0.1J, 
We are also given that 


J, =107210+0.222J, 


0.1J, = 107214 0.0222, 


Hence 
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J, = 42.5x10° x 6004 + 71035 + 10721 + 0.0222/, 
0.97778, =5508+81756 

J, =89247 W/m? 

Also 

E,g—J, _ 56.7 10° x 600* — 89247 


ae 1/3 
Ey 


= —246x10° W/m? 


qy = 


Negative sign indicates J, >F', > E, , so net flux is into the piston crown. 


c) To make the analysis more realistic, it needs to be extended by including convection from the 
piston crown, and cylinder liner. Radiation from the piston underside also needs to be included. 
We then carry out analysis over a complete engine cycle. 


Example 4.5 


The figure below shows the variation of view factor F;,; with geometric parametres h / L and W/ L for 
the case of two rectangular surfaces at right angles to each other. This plot is to be used to model the 
radiative heat transfer between a turbocharger housing and the casing of an engine management 
system. The horizontal rectangle, W = 0.12 m and L= 0.2 m, is the engine management system and is 
denoted Surface 1. The vertical rectangle, h = 0.2 m and L = 0.2 m, is the turbocharger casing and 
denoted by Surface 2. The surrounds, which may be approximated as a black body, have a temperature 
of 60°C. 


a) Using the graph and also view factor algebra, evaluate the view factors: F) 2, Foi, Fi; and F.3 


b) By applying a grey-body radiation analysis to Surface 1 with ¢; = 0.5, show that the radiosity 
Jy is: 


J, = 28.35 x 10° T,* + 0.135 Jy + 254 (W/m’) 
The following equation may be used without proof: 


E,,-J,  & 
Weave I) 
i i J= 
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c) A similar analysis is applied to Surface 2 with € = 0.4 obtained the result: 
Jo = 22.7 x 10° T>* + 0.097 J, + 350 (W/m’). 


Use this to estimate the surface temperature of the engine management system when the turbocharger 
housing has a surface temperature of T, = 700K. 


W /L=0.05 


0.1 04 0608 | 2 4 10 
h/L 
Solution 
12 

Oe. PRUE 0.2m 
L 0.2 L ie 
From the figure: F,, = 0.27 0.12 2m 
AF, = A), 

A 2 
| Neeley nee yt TS) 

A, h 0.2 


Fy t+Fy +h3 =1 


FF, =0 
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F,, =1-F,, =1-0.27=0.73 


FP, +f, + Fy, =1 


F,, =0 


F,, =1-F,, =1-0.162 = 0.838 


For a grey body radiative heat transfer in an enclosure (n surfaces) 


1 
og, MIF 


E. 


L 


Applying for surface 1, 1= 1 (the casing) 
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E,,-J 
Fil —J,)+ Fy (J, -J;) 
1 
é) 
£,, =oT, 
Joel, 
l-e, 7 1-0.5 -10 
&; 0.5 
So 


= ol; + FJ, +F,oT, 


J, 
I+F, +f; 
7p 287% 10° 7? +0.27J, +0.73x 56.7 x10” x 3334 
1+0.27+0.73 
J, =28.35x10°T? +0.135 J, +254 W/m? (1) 
c) 


Given: J, =22.68x10°T7, +0.0972 J, +350 W/m? 

J, =22.68x10° x 7007 + 0.0972 J, +350 W/m? 

J, =5796 + 0.0972 J, (2) 
Substituting from equation 2 into equation 1: 

J, = 28.35x10°T* +.0.135 (5796 + 0.0972 J,)+254 W/m? 

Which gives: 

J, =28.7x10°7,7 +1050 W/m? 

Applying a heat balance to surface 1 


qin = out 
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Ee md 
G_, =—-| 2 + | =-57.9x 10° T — 28.7x10°T* +1050 
le; 


é; 


Yin = —28.x10°T,* +1050 
dow = €,0(T* -T3)=0.5x56.7x10 (71 - 333") 


Combining and solving for T), gives: 


T, =396 K 


Note that gin =-q since q is out of the surface wheng>0 . 
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5. Heat Exchangers 


Example 5.1 


A heat exchanger consists of numerous rectangular channels, each 18 mm wide and 2.25 mm high. In 
an adjacent pair of channels, there are two streams: water k = 0.625 W/m K and air k = 0.0371 W/m 
K, separated by a 18 mm wide and 0.5 mm thick stainless steel plate of k = 16 W/m K. The fouling 
resistances for air and water are 2 x 10* m’ K/W and 5 x 10° m’ K/W, respectively, and the Nusselt 
number given by Nupn= 5.95 where the subscript 'D,' refers to the hydraulic diameter. 


a) Calculate the overall heat transfer coefficient ignoring both the thermal resistance of the 
separating wall and the two fouling resistances. 

b) Calculate the overall heat transfer coefficient with these resistances. 

c) Which is the controlling heat transfer coefficient? 


Solution: 


Hydraulic Diameter = 4 x Area / Wetted perimetre 


3 3 
D, = 4x 225%10 «18x10 ai 
(2.25 +18) x 107 


ju Nu yk 
D, 
(@) agg = 08” _ O30 1m? K 
4x10 
air eI! = 55.186W /m?K 
4x10 


-1 
U= = : =52.1 W/mK 
930 55.186 


-1 
b) y =|(Sx Legh sa" =50.2 W/mK 


c) The controlling heat transfer coefficient is the air heat transfer coefficient. 
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Example 5.2 


A heat exchanger tube of D = 20 mm diameter conveys 0.0983 kg/s of water (Pr = 4.3, k = 0.632 W/m 
K, p = 1000 kg/m’, 1: = 0.651 x 10° kg/ms) on the inside which is used to cool a stream of air on the 
outside where the external heat transfer coefficient has a value of h, = 100 W/m’ K. Ignoring the 
thermal resistance of the tube walls, evaluate the overall heat transfer coefficient, U, assuming that the 
internal heat transfer coefficient is given by the Dittus-Boelter relation for fully developed turbulent 
pipe flow: 


Nu, = 0.023Re{* Pr’ 


Solution: 


_ pVD— 4ni 4x 0.0983 


- - = 9613 
uw = mDu 2x0.02x0.651x107 


Nu, = 0.023 x 9613°* x 4.3°% = 63 


hD 
ers 
Fos Nu pk _ 63.3x 0.632 “ANNU 
D 0.02 
1 17° 
U= aortas | =95.2W /m°K 
2000 100 
Example 5.3 


a) Show that the overall heat transfer coefficient for a concentric tube heat exchanger is given by the 
relation: 
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With the terminology given by the figure below 


b) A heat exchanger made of two concentric tubes is used to cool engine oil for a diesel engine. The 
inner tube is made of 3mm wall thickness of stainless steel with conductivity k = 16 W/m K . The 
inner tube radius is 25mm and has a water flow rate of 0.25 kg/s. The outer tube has a diameter of 
90mm and has an oil flow rate of 0.12 kg/s. Given the following properties for oil and water: 


oil: 
C, =2131 JkgK, m=3.25x10~ kg/ms, &=0.138 W/mK 
Water: 

C, =4178 JkgK, =725x10° kg/ims, k=0.625 W/mK 
Using the relations: 


Nu, =5.6 Re, < 2300 
Nu, = 0.023 Re?* Pr°* Re, > 2300 
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Heat Exchangers 


Calculate the overall heat transfer coefficient. 


Which is the controlling heat transfer coefficient? 


If the heat exchanger is used to cool oil from 90°C to 55°C, using water at 10°C calculate the length of 


the tube for a parallel flow heat exchanger 


Solution: 

a) 

For the convection inside 

O= AA(T,-T)) 

OQ =2nr, Lh(T,-T)) 

For the convection outside 

OS Ah, 1) 

OQ =2nr, Lh,(T,-T,) 

For conduction through the pipe material 


aT 
=—22 rk — 
g " dr 


(1) 


(2) 
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(oft 
2QarL)r 


Integrating between | and 2: 


r,-T,=|—£ 
2ar, Lh, 


Adding 4, 5 and 6 


7-2 [etn 11 ) 


“Ot| i he ae 


i 


Rearranging 


b) 


(3) 


(4) 


(5) 


(6) 


i) To calculate the overall heat transfer coefficient, we need to evaluate the convection heat transfer 


coefficient both inside and outside. 
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Re= pV ,,D 
L 
For water: 
. 2 
Vn = oe A = a 
fA 4 
Ree 4m _ 4x0.25 : — 978] 
my 2mx0.05x725x10 
C -6 
Pr= HC, Zz 725x10~ x 4178 ~ 485 


k 


0.625 


Re > 2300 (turbulent flow) 


Therefore: Nu, = 0.023Re‘s* Pr°* = 0.023 x8781"* x 4.85°% = 62 


What do you see? 


A skateboarder? 


We see an All Options Junior Trader leaving 
the office after a successful day. 
He has the freedom to make his own decisions. 
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Nuk — 62x 0.625 
D 0.05 


=775 Wim K 


From which: h, = 


For oil: 


D= 4 Area _ An(r; -1r7) 


[_=———_— == = 2(r, — 7) = 2(0.045 — 0.025) = 0.034 m 
Perimeter 2z(r, +1r,) 


pV,,D, _ 2m(%, -1,) 2m 2x0.12 
Re — m 1 — a — — 
L mn—re)u ar, +r, 2x(0.045+0.028)x3.25x107 


= 33 


Re < 2300 (Laminar flow) 
Therefore: Nu, =5.6 


— Nuyk _ 5.6x0.138 


h, = 
D, 0.034 


=22,1 W/m2 K 


-1 
Ole le | ae. Wark 
16 (25) 725x0.025 22.7 


ii) The controlling heat transfer coefficient is that for oil, A, because it is the lower one. Changes in 
h, will cause similar changes in the overall heat transfer coefficient while changes in h; will cause 
little changes. You can check that by doubling one of them at a time and keep the other fixed and 
check the effect on the overall heat transfer coefficient. 


iii) 7,,=90°C, 7,,=10°C, 7,,=55°C 
T., is unknown. This can be computed from an energy balance 


For the oil side: 


QO = 1i1,C,,, (Ty, —T,,) = 0.12 x 2131(90—35) = 8950 W 
OQ =1h,C,,(L,, -T,) = 0.25x4178(T,, —10) = 8950 W 


Therefore T,, =18.56°C 


Evaluate LMTD 
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AT, = 90-10 = 80°C 
AT, = 55—-18.56 = 36.44°C 


AT = Se ee eee 
In(AT,/AT,)  1n(36.44/80) 


QO =UAAT 


Im — U, ‘“ 2rr, LAT,, 
= g = one =41.5m 
U,x2ar,AT,, 21.84x 27x 0.028 x 56.1 


m 


Example 5.4 


Figure (a) below shows a cross-sectional view through part of a heat exchanger where cold air is 
heated by hot exhaust gases. Figure (b) shows a schematic view of the complete heat exchanger which 
has a total of 50 channels for the hot exhaust gas and 50 channels for the cold air. The width of the 
heat exchanger is 0.3m 


Using the information tabulated below, together with the appropriate heat transfer correlations, 
determine: 


i. the hydraulic diameter for each passage; 
li. the appropriate Reynolds number; 
iil. the overall heat transfer coefficient; 
iv. _ the outlet temperature of the cold air; 

v. and the length L 


Use the following relations: 


Using the relations: 


Nu, =4.6 Re, < 2300 
Nu, = 0.023Re?* Pr‘? Re, > 2300 
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Data for example 4.4 


H=3mm 


Cold air > 
Hot gas <—______ 3mm 


Cold air —________, 

Hot gas ¢—$£§.- 

Cold air —_______, t 
Hot gas t 


(a) Cross sectional view through part of the Heat Exchanger 


W=0.3m 


(b) Schematic diagram of the complete heat exchanger 


Hot exhaust inlet temperature 100°C 
Hot exhaust outlet temperature 70°C 
Cold air inlet temperature 30°C 
Hot exhaust total mass flow 0.1 kg/s 
Cold air total mass flow 0.1 kg/s 
Density for exhaust and cold air 1 kg/m® 
Dynamic viscosity, exhaust and cold air 1.8x10-5 kg/m s 
Thermal conductivity, exhaust and cold air 0.02 W/m K 
Specific heat capacity, exhaust and cold air 1 kJ/kg K 
Heat exchanger wall thickness 0.5mm 
Heat Exchanger wall thermal conductivity 180 W/m K 
Hot exhaust side fouling resistance 0.01 K m7/W 
Cold air side fouling resistance 0.002 K m7/W 
Solution: 
Ree 

L 


L=D, — (Hydraulic diameter) 


D 


_ 4xcrosssectionalarea 4xwxH — 4x0.003x0.3 


a = =5.94 mm 


h 


perimenter 


2(w+H)  2(0.003 + 0.3) 
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For a single passage: 


(a/50) ___0.1/50 
(H x w)p (0.003 x 0.3) 1 


Y= =2,22 mis 


_1x2.22x5.94x107 


Re 
1.8x10> 


= 733 


Re < 2300 (laminar flow) 
Nu, = 4.6 


— Nupk — 4.6x0.02 


b. See 1g 715.5 Wi m'K 
h . x 


h 
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Since the thermal properties are the same and the mass flow rate is the same then the hot stream and 
cold stream heat transfer coefficients are also the same. 


Note that if the third term in the brackets that includes the resistance through the metal is neglected, it 
will not affect the overall heat transfer coefficient because of the relatively very small thermal 
resistance. 


Q = mC (7, =1 5.) = mC, (7. =) 


Tas =T,,; +(7,,; —T,.,) = 30+ (100-70) = 60°C 


Also 
Q = UAAT,, 
T,,, 1S constant in a balanced flow heat exchanger 


Tm = 100 — 60 = 70 —30 = 40°C 


Im 


Q=mC, (7, = ) = ~* 1000(100-70)=60 w/passage 


Area of passage: 
4-2 = am? 
UAT,, 7.1x 40 


And since: 4=wxZ 


_ 0.211 
0.3 


L = 0.704 m 
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